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HEAT CAPACITY LAG OF GASEOUS MIXTURES 

By Thomas D. Rossing, Robert C. Amme, 
and Sam Legvold 


SUMMARY 


The relaxation times for the excitation of molecular vibrations in 
a number of heavy gases have been measured with an acoustic interferome- 
ter. All the gases studied were found to have a single relaxation time, 
indicating that intermodal coupling in these gases is strong. 

Only binary collisions are found to be important in the excitation 
of vibrations. The probability, in a given collision, of exciting or 
deexciting molecular vibrations in the halogen- substituted methanes 
appears to depend upon the relative energy of approach of the colliding 
molecules rather than upon the relative velocity. 

The excitation of vibrations in polyatomic gases may well involve 
the formation of a complex molecule with a very short life. 


INTRODUCTION 


About 30 years ago, investigators in the field of ultrasonics dis- 
covered that the velocity of sound in polyatomic gases is a function of 
the sonic frequency. This phenomenon is known as the dispersion of 
sound and may be explained by the failure of the translational, rota- 
tional, and vibrational energies of gas molecules to maintain thermal 
equilibrium with each other at high acoustic frequencies. An irregularly 
high acoustic absorption accompanies this dispersion in acoustic velocity. 

Numerous investigators have measured the dispersion of sound caused 
by the lagging of the vibrational energy (also referred to as the thermal 
relaxation process). Nearly all the investigations have been made in 
diatomic and triatomic gases in which the vibrational relaxation time is 
relatively long and the region of dispersion is at easily attainable 
frequencies. Advances in electronics have recently made possible the 
accurate measurement of acoustic velocities at very high frequencies. 
Considerable interest has been aroused in the study of dispersion due to 
the lag of rotational and translational energies; study of the vibrational 
energy lag in heavy organic gases, which also necessitates high-frequency 
measurements, has been neglected. It is the purpose of this paper to 



2 


MCA TN 3558 


report some measurements of sonic dispersion made in heavy gases and to 
attempt clarification of certain fundamental properties of the vibrational 
excitation process in these gases. 

This work was conducted at Iowa State College under the sponsorship 
nnri with the f inan cial assistance of the National Advisory Committee for 
Aeronautics . 

SYMBOLS 



^atomic 

Co 

Crot 

Cfrans 

C v 

C vib 

C (o 


C 


00 



h 

k 


Beattie-Bridgman gas constants 

contribution to Cy made by atomic excitation 

molar heat capacity at constant volume at frequencies 
below dispersive region 

contribution to C v made by molecular rotations 

contribution to Cy made by molecular translational 
motion 

molar heat capacity at constant volume 

contribution to C v made by molecular vibrations 

effective Cy at angular frequency cd in dispersive 
region 

molar heat capacity at constant volume at frequencies 
above dispersive region 

velocity of light 

sonic frequency 

sonic frequency at which inflection point of dispersion 
curve occurs 

Planck’ s constant 

Boltzmann's constant 
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M 

N 


n 


n 0 

n l 



R 

s 

T 

t 

V 

V bb 

V ideal 

Vo 

V M 


v 


Z 


10 


molecular weight 

average number of collisions experienced by a molecule 
per second 

ratio of e* to hv-^ 

number of molecules of gas per unit volume 

number of molecules per unit volume whose vibrations are 
in ground state 

number of molecules per unit volume whose vibrations are 
in first excited state 

average of p^^ over proper velocity or energy 
distribution 

pressure 

probability that a molecule- in first excited vibrational 
state will be deexcited by a collision; a function of 
velocity or energy of approach 

universal gas constant 

range of intermolecular forces 

temperature, °K 

time, sec 

velocity of sound 

velocity of sound in a Beattie-Bridgman gas 
velocity of sound in an ideal gas 


velocity of 


sound at frequencies below dispersive region 


velocity of sound at frequencies above dispersive region 


relative velocity of two colliding molecules 

average collision lifetime of first excited state; aver- 
age number of collisions necessary to deexcite a mole- 
cule in first excited state 
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e 

e* 

0 

P 

v 

v 

V 1 

I 

P 

P 

a 

aj 


“inf = 2rtf inf 


steric factor; probability that a collision ■will take 
place with molecules favorably oriented for excitation 

relative approach energy of two molecules 

mi ni mum value of e required for vibrational excitation 

vibrational relaxation time 

reduced molecular weight of two colliding molecules 
vibrational frequency per second 
wave number per centimeter, • v/c 
vibrational frequency per second of lowest mode 

velocity correction due to absorption 

density, g/liter 

density, moles /liter, p/M 

see equation (27) 

angular frequency of sound, 2rtf 


THEORY OF ACOUSTIC DISPERSION 
Propagation of Sound 


Richards (ref. 1, 1939) bas derived an expression for the velocity 
of sound in an isotropic homogeneous medium: 

= +(|gf “ (1) 

Vdp/T VST/p 

in which p is the pressure; T, the temperature; Cy, the molar heat 
capacity at constant volume; p, the density; and M, the molecular 
weight. For the velocity in an ideal gas, the derivatives may be eval- 
uated from the ideal-gas equation p = ^ = ^=jp. and substituted into 
equation (l): 
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'ideal 



( 2 ) 


An equation of state frequently applied to polyatomic molecules is 
that of Beattie and Bridgman: 


p - - 4 (3) 

where 

A = (l - |)a 0 
B = (l - £>0 

a, b, A Q , and B 0 being constants characteristic of the given gas. 
This equation may also be written in terms of p = 

p = ETp + (RTB 0 - Aojp 2 + ^A Q a - RTB 0 bjp 5 + . . . (4) 

Substitution of equation (k) into equation (l) produces an expression 
for the square of the velocity in a Beattie -Bridgman gas : 

V bb 2 - sj™ + 25(^0 - A o) + 55 2 ( A o a - KBVO] + 

+ 2pB 0 + 3 2 (b o 2 - 2B 0 b)] (5) 

Terms in p3 and higher have been omitted. Note that as p becomes 
very small (low pressures) equation (5) reduces to equation (2). 

In any gas molecule, the specific heat at constant volume may be 
expressed in terms representing contributions from translatory motion, 
internal vibrations and rotations, and atomic excitation: 


— ^trans + ^rib + ^rot + ^ atomic 


(6) 
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The value of each quantity depends upon temperature, and, according to 
the equipartition theorem of energy, in a body at equilibrium the energy 
(and therefore the heat capacity) must be distributed among these terms 
in a definite manner. At room temperatures most atoms are in their 
ground state and C a . hrmn . p is negligible. If the equilibrium between 

these types of internal energy is disturbed, each mode again approaches 
its equilibrium value at a rate characterized by its own relaxation time. 

Vibrational relaxation times cover the range 10“^ to 10 second. 

As long as thermal disturbances occur slowly, C is correctly given 
by equation (6). Since the period of a disturbance is so shortened as to 
be of the order of one of these relaxation times (as happens when a sound 
wave of high frequency is incident), the equilibrium process is unable to 
follow the rapid thermal cycle, and the net effect is that a part of the 

heat capacity drops out completely. Thus Cy and depend upon fre- 

quency throughout the dispersive region. 

If C Q denotes the heat capacity at constant volume below the dis- 
persive region and C^, the heat capacity above this region, then 


C = 

CO 


■'trails 


+ C. 


rot 


( 7 ) 


and 


C o = C oo + °vib (8) 

Following the technique of Richards (ref. 1, 1939), the heat capacity at 
acoustic angular frequency o> within the dispersive region may be written 


Coi - + 


C o - Cco 

1 + io>e 


(9) 


where 0 Is the appropriate relaxation time and cd is 2n times the 
acoustic frequency f. Substitution of equations (7) and (8) into equa- 
tion (2) yields 





( 10 ) 


( 11 ) 
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Similarly, substitution of equation (9) into equation (2) will yield a 
complex expression, the real part of which represents the phase velocity 
of sound at frequencies within the dispersive region: 


V. 


ideal 


RT 

M 


1 + R 


C Q + a&gfiX 

C 0 2 + o&Zcj) 


( 12 ) 


The inflection point of equation (12) may be found by differentiation: 


cjd. 


'inf 


ec 


(13) 


Similarly, substitution of equation (9) into equation (5) gives an 
expression for the square of the velocity in a Beattie-Bridgman gas: 



The dispersion curve for an ideal gas, plotted from equation (12), is 
shown in figure 1. 


Multiple Dispersion 

Thus far it has been assumed that for a given gas the relaxation 
process may be characterized by a single relaxation time 9 . This is 
equivalent to assuming that all modes of vibration are excited with the 
same relaxation time and, furthermore, that al 1 vibrational transitions 
occur with the same relaxation time as transitions from the ground state 
to the first excited state. That all transitions between the states of 
a given mode of vibration have the same relaxation time has been shown 
by Landau and Teller (ref. 2, 1936). If, however, the couplings between 
different modes of vibration were so weak- that they must be excited inde- 
pendently by collisions and therefore possess different relaxation t ime s, 
a step-shaped velocity-dispersion curve would result. 

Some evidence favoring the existence of multiple dispersive regions 
has been reported. Pielemeier (ref. 3, 1943) and Busc hmann and Schafer 
(ref. 4, 1941) have reported a possible second step on the velocity- 
dispersion curve of carbon dioxide. Alexander and Lambert (ref. 5, 1942) 
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found three dispersive regions in acetaldehyde. A later article by 
Lambert and Rowlinson (ref. 6, 1950) , however, points out that the mul- 
tiple dispersions of velocity are due to errors in calculation resulting 
from an ungrounded assumption that dimerization equilibrium does not keep 
pace with the acoustic cycle. 

Great difficulty would be experienced in the observation of multiple 
dispersive regions in triatamic molecules, if indeed such exist. At room 
temperature one mode usually accounts for most of the vibrational contri- 
bution to specific heat. Any steps on the velocity-dispersion curve due 
to the other modes would be so small that they would be difficult to 
detect. Detection of multiple dispersion due to a lack of intermodal 
coupling, if it exists, might be most probable in molecules such as the 
halogen-substituted methanes since these have modes of low frequency 
which can be excited at room temperatures . On the other hand, if inter- 
modal coupling is strong enough so that energy is transferred readily 
from one mode to another, only a single relaxation time will be observed. 
Presumably this will be the relaxation time for the mode of lowest 
frequency. 


Trace Catalyst Effect 

A phenomenon which has been observed in relaxation-time studies is 
the "trace catalyst" effect. Collisions between unlike molecules are 
sometimes much more efficient in the excitation of vibrations than those 
between like molecules. Eucken and Becker (ref. 7, 1934) observed that 
in a collision between two molecules which have a chemical affinity for 
each other the probability of vibrational excitation is large. Landau 
and Teller (ref. 2, 1938) suggested that excitation is more probable 
when the two colliding molecules approach each other at a high velocity, 
thus explaining the effectiveness of hydrogen and helium as catalysts. 
Several, important catalytic gases are discussed by Walker (ref. 8, 1951)* 


Multiple Collisions 

In general, a certain average number of collisions are required to 
excite or deexcite the molecular vibrations in a given gas. Since the 
number of binary collisions is proportional to the pressure of the gas, 
it follows that the relaxation time will be inversely proportional to 
the pressure when binary collisions are responsible for vibrational 
excitation. If triple collisions are primarily responsible, the relaxa- 
tion time will be inversely proportional to the square of the pressure. 

In equations such as expressions (12) and (l4) dealing with acoustic 
dispersion, the independent variable is always c 00 rather than cd. 

Hence, if the relaxation time 0 is inversely proportional to the 
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pressure, the dispersion curve may well "be drawn with log f/p as the 
abscissa instead of log f . This is a great advantage in experimental 
work, since it provides a method of obtaining a continuous range of fre- 
quencies without changing the frequency of the sound source. Of course, 
it must first be ascertained which type of co lli si on is responsible for 
vibrational excitation so that the correct pressure dependence is used. 

No conclusive evidence has been presented supporting any importance 
of collisions of orders above the second. Walker (ref. 8, 1951) suggested 
that certain data of previous investigators, involving carbon oxysulfide 
with added quantities of helium and argon, might be interpreted on the 
basis that triple collisions are of greater importance than the binary 
type. Walker, Rossing, and Legvold (ref. 9* 195^) later reported that 
triple coll i sions are not important in the excitation by helium, argon, 
and nitrogen of nitrous oxide, which is very similar in structure to 
carbon oxysulfide. 


THEORY OF VIBRATIONAL EXCITATION BY COLLISIONS 
Relaxation Theory 


A gas will be considered which has n^ molecules with vibrations 
in the first excited state and n^ molecules in the ground state per 
unit volume . If at all times n Q - n-j_ = n, a constant, the continuous 

exchange of molecules between states 0 and 1 may be characterized by the 
equation 

dnn 

= f oi n o " f io n i ^ 

in which f^ specifies the number of 0 — >1 transitions per ground 
state molecule per second and f^Q, the corresponding number of 1 — >0 
transitions per excited molecule. 

dn-i 

When the gas is in equilibrium, — — = 0 and, therefore, according 

dt 

to the principle of detailed balancing, the total number of transitions 
in each direction are equal. At equilibrium 

-hv/kT 

“l = ^ 



and so 
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< 

-hv/kT 

f 01 = f io e 


(17) 


Tn "the above equations li and k are the Planck and Bolt zm a nn constants , 
respectively, v is the frequency of molecular vibrations, and T ip the 
temperature in °K. 


Since n = n^ + np at al 1 times, equation (15) may be written: 


dn-3_ 


dt 


+ 



(18) 


Since n is constant, the solution to this differential equation for 
nonequilibrium is 


n^ = c^ + C2e 


"( f oi +f io) t 


(19) 


where c^ and are constants. The relaxation time of this process 

can then be seen to be 


0 


1 

f 01 + f 10 


( 20 ) 


Here 0 is the time necessary for the departure from thermal equilibrium 
to be reduced to l/e of its initial value . Some authors assume that 
hv/kT » 1 and f Q1 « f 10 and, hence, that 9 » 1 jf-^Q . At vibrational 


frequencies less than about v = 1.5 


13 -1 

x 10 second 


(corresponding to 


v = 500 centimeter ^ in spectroscopic wave numbers, since v = v/c), 
this assumption is ungrounded. 


If the average number of collisions which a molecule ma k es per sec- 
ond is denoted by N, f Q1 and f^o may be written in terms of Pq^ 

and P 1Q , the probabilities that a molecule in the ground or first excited 
state will be excited or deexcited, respectively, in a collision. Then, 


f 01 « P 0lN 

f 10 = P 10N 


( 21 ) 
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The reciprocals of P Q1 and P 10 , Zqi and will give the mean 

collision lifetimes of the ground and the first excited states, respec- 
tively. Substitution of equations (17) and (21) into equation (20) 
produces 


0 


NP 10 (l + e " hv / kT ) 


1 

BP 0 l(l + e 11 ^) 


Therefore : 


or 


J 0i 1 + e 


hv/kT 


Z 10 1 + e -^v/M 


Z 01 “ Z 10 e 


hv/kT 


In case hv « kT (not realizable at ordinary temperatures) 


( 22 ) 


( 23 ) 


( 2*0 


Z 01 *** Z 10 ( 25 ) 


Prediction of P 

10 

Landau and Teller (ref. 2, 1936) have developed a theory which has 
been very successful in explaining energy-transfer probabilities in 
relatively simple molecules. They write an expression for the probability 
of deexcitation of a vibrating molecule p^^ in a single collision when 

the relative velocity of approach of the molecules lies within the 
range v to v + dv: 


P 


10 


-2rtvs/v 

ae 


(26) 


where a is a steric factor which gives the probability that a collision 
will find the molecules favorably oriented for excitation or deexcitation. 
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s is the range of intermole cular forces, and v, the vibrational fre- 
quency. Bethe and Teller (ref. 10, 19^-2) integrate the above expression 
over a Maxwellian distribution of velocities to obtain 



(27) 


where 


a 


0 . 090 (vs ) 


-_x 2 /3/u \l/3 


¥ 



p is the reduced molecular weight of the colliding molecules, v is 
the vibrational frequency in centimeter - -*-, measured in units of 

10 - 9 centimeter, and T is in °K. 

Recent theory by Schwartz, Slavsky, and Herzfeld (ref. 11, 1952), 
based on the approach of Zener, makes use of an exponential repulsion 
and a one -dimensional model. These authors arrive at a somewhat more 
complicated expression for P^q which, they point out, has essentially 

the same dependence upon p, v, and T as the theory of Landau 
and Teller. 

Fogg, Hanks, and Lambert (ref. 12, 1955) have measured the vibra- 
tional relaxation times of several halomethane vapors and have calcu- 
lated the number of collisions necessary to excite the lowest modes of 
vibration. At least four of these gases appear to fit an empirical rela- 
tion of the form Pq-^ = e in which v-^ is the frequency of the low- 

est vibrational mode and X is a proportionality constant characteristic 
of the particular group of gases. This suggests that the probability of 
vibrational excitation depends upon the energy of approach of the col- 
liding molecules rather than upon the velocity of approach. An energy 
excitation theory, based on the above suggestion, will now be furt he r 
developed by examination of the deactivation process. 

It may be assumed that the probability ' p 1Q of deexciting a vibra- 
ting molecule in a given collision can be written in terms of the rela- 
tive energy of approach e of the colliding molecules: 
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p = a 

10 

All 

£* 

= 0 

6 < 

e* 


(28) 


In equation (28), a is the steric factor arid g* is some minimum 
relative energy of approach which, must he exceeded in order to produce 
deexcitation. If it is assumed that the mode of lowest frequency 

is the first to he deexcited, e* may he written e* = nhvj_ where n 
is some constant. 


Integration of equation (28) from e = g* to e = <» gives the 
prohahility of deexcitation in a single collision: 




= ae 


-nh\q_/kT/^ lv i \ 
\ kT + 7 


(29) 


The corresponding prohahility for excitation may he written: 


p 01 . + x ) ( 50 ) 

The above relations picture vibrational excitation and deexcitation 
as involved processes requiring a certain minimum energy of approach. 

It is tempting to assume that vibrational excitation and deexcitation 
take place only when a molecule is able to penetrate very deeply into 
the force field of another, the penetration requiring a certain amount 
of energy . It is even possible that the collision partners form a com- 
pound molecule with a very short life, during which energy may pass 
readily between translational and vibrational degrees of freedom. The 
process of formation of the compound molecule would require a certain 
minimum amount of energy. It should be noted that the probabilities 
for the excitation and deexcitation processes differ by a factor 

e-hVi/kT. Fogg, Hanks, and Lambert (ref. 12, 1953) have apparently 
treated them as being the same . 


Trace catalysts may still exist for gases which obey the energy 
excitation theory, but gases which catalyze diatomic and triatomic gases 
may not produce the same result in polyatomic gases. Hydrogen and 
helium apparently are strong catalysts in diatomic and triatomic gases 
because they have a very high velocity at ordinary temperatures. To be 
a catalyst in a gas obeying the energy excitation theory, however, a 
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gas must tie able to penetrate the force field of another easily, or to 
form a compound molecule readily, and have a fairly large steric factor. 


Heat Capacities of Gases 

According to equation (6), the heat capacity of a gas is made up of 
contributions from molecular vibrations, rotations, and translations. 

At sufficiently high temperatures, each degree of translation or rota- 
tion contributes ^R, and each vibrational mode contributes R times the 
degree of degeneracy. 


At room temperatures, C ro ^. is fully developed in polyatomic gases. 
Thus the rotational contribution from linear polyatomic molecules is R, 
anil the contribution from nonlinear polyatomic molecules is ^R, since 

they have moments of inertia along three perpendicular directions. 
Translational energy, also characterized by three degrees of freedom, 

contributes ?R. 

2 


Vibrational heat capacities, on the other hand, are rarely complete. 
Assuming that the molecules are perfect harmonic oscillators, the vibra- 
tional contribution to heat capacity may be written in terms of the fre- 
quencies of the n modes each having a degeneracy of order j: 


n 

Wib ~ S - 
i=l 



(3D 


Spectroscopic data indicate that molecular vibrations are not per- 
fectly harmonic, and thus a transition from the second excited state to 
the first excited state does not produce the same frequency as that pro- 
duced by a transition from the first state to the ground state. The 
error in V G introduced by the neglect of this anharmonicity in calcu- 
lating is about 0.03 percent. The total heat capacity at constant 

volume, at room temperature and below the dispersive region in polyatomic 
molecules, may then be given as 


C 0 = 3R + (32) 

where C^b is calculated from equation (31)- Above the dispersive 
region 


C = 3R 
00 


( 33 ) 
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EXPERIMENT 

Apparatus 


The velocity of sound has been measured by means of an interferometer 
of the Pierce type . A cross section of the acoustic chamber used to meas- 
ure wavelengths is shown in figure 2. A gold-plated X-cut quartz crystal 
is clamped in its nodal plane and is driven near its frequency of reso- 
nance. An optically flat stainless- steel reflector is connected by a 
precision ground invar rod to a micrometer screw which has a least count 
of 0.001 millimeter. The crystal and the reflector are maintained paral- 
lel within 0.1 wavelength by means of three leveling screws. 

The acoustic chamber is constructed from stainless steel, nl 1 seams 
being silver soldered. A copper gasket insures a vacuumtight seal. Thp 
invar rod, a mercury thermometer, and the osc illat or leads pass through 
0-ring vacuum seals. The entire chamber is immersed in a constant- 
temperature water bath. 

The crystal oscillator circuit is shown in figure 3- It is a varia- 
tion of the Pierce type of circuit, using precision components throughout, 
and drives the crystal at a very constant frequency and amplitude. Thp 
theory of its design is discussed by Rossing (ref. 13, 1952). The fre- 
quency of the crystal is determined to a precision of 1 part in 10,000 
by heterodyning the output of the oscillator with the signal from some 
convenient broadcast station to produce a beat note in the audiofrequency 
range. Most of the measurements were made with a 300-kilocycle crystal, 
although a 1,000-kilocycle crystal was sometimes employed. 

Sharp peaks in the input impedance of the crystal indicate the 
positions of the reflector for which standing waves result. The cir-. 
cult used to measure small changes in the crystal impedance is shown in 
figure 4. Voltages from the crystal are isolated by two cathode fol- 
lowers, Independently rectified and fed into an unbalanced direct-current 
amplifier. A biased microammeter with five ranges measures the changes 
in the output of the direct-current amplifier. A General Electric vacuum 
thermocouple, which generates a voltage proportional to the amount of 
radiofrequency current through it, is used to measure the crystal cur- 
rent. The impedance -measuring circuit and the crystal oscillator receive 
their power from a stable power supply with an electronic voltage regu- 
lator of the degenerative feedback type. The power-supply circuit is 
shown in figure 5> 

The gas -handling system is shown schematically in figure 6. A two- 
stage mechanical vacuum pump is capable of evacuating the entire system 
to a pressure of 0.2 micron. The pressure in the acoustic c hamb er is 
measured by a two-meter manometer and McLeod vacuum gage. A drying 
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chairiber is available to remove traces of water- vapor impurity from gases; 
test gases can be drawn into a liquid-air cold trap and distilled into 
the c hamb er to remove traces of noncondensahle gases. Kerotest packless 
diaphragm valves are used throughout the system. Before each measurement, 
the entire system is evaluated and out gassed to assure a gas sample of as 
high purity as possible. 


Experimental. Procedure 

The velocity of sound was measured in each gas at a variety of 
f/p ratios by means of the acoustic interferometer. The corresponding 

values of V 2 were then plotted on semilogarithmic graphs in figures 7 
to 20; v ia ea q 2 was plotted using equation (12) with an arbitrary value 

of 0. This curve was then shifted along the f/p-axis until it best fit 
the experimental data, since according to equation (12) the effect of 
changing 0 is merely to shift the entire curve along this axis. In 
those gases for which the dispersive region occurred at low frequencies 
and for which the Beattie-Bridgman constants were known, the V-^ 2 curve 

has been drawn using equation (l^) . The experimental values of V 2 were 

2 2 

then corrected by adding the value of - V-^ measured along the 

f/p ordinate corresponding to that point. Nearly all the values of V 2 
axe the result of 10 independent determinations of wavelength. 

Since C rot is fully developed in polyatomic gases at room tempera- 
tures, C m was calculated from equation (33); C vib was tlien calculated 

using the vibrational wave numbers observed spectroscopically and equa- 
tion (32); C Q was then computed by equation (32). Thus the anharmonicity 
of the molecular vibrations was taken into account. Values of are 

given in table I; also given are values of C-^, that part of con- 

tributed by the lowest vibrational mode. 

The gases used in the experiment, which were of the highest purity 
available, were as follows: CHCTF^ (chlorodif luoromethane), CHCl^F 

(dichlorofluorometbane), CCI 2 F 2 (ddchlorodif luoromethane ) , CH^Cl 
(chlorome thane), CH^Br (br omome thane ) , CCljF (trichlorofluoromethane), 

and CCIF^ (chlorotr if luoromethane), which were purchased from Matheson Co., 
Inc., ana CBrF^ (bromotrif luoromethane ) , CH^CIF ( chlorof luoromethane ) , 

CH 2 F 2 (dif luoromethane ), CBr 2 F 2 (dibromod if luoromethane ) , CBrCU?^ 
(bromochlorodif luoromethane ) , CFi,. (tetraf luoromethane ) , and CHF^ 
(trifluorometbane), which were furnished by the Jackson Laboratory of 
E. I. du Pont de Nemours & Co., Inc. The purities reported by the sup- 
pliers of the gases are listed in table I. Several of the above gases 
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were stored in cylinders as gas over liquid. These gases were withdrawn 
in their liquid plaase to minimize contamination 'by air anri noncondensable 
gases. Others were liquified by use of the liquidinitrogen cold trap and 
then evaporated into the system after pumping off al 1 noncondensable 
vapor. The close fit of the low-frequency exper ime ntal data to the cal - 

p 

culated values of V 0 demonstrates the high degree of purity of the 
gases used. 


Values of were taken from the curves in figures 7 to 20. The 

corresponding values of 9 were then calculated by use of equation ( 13 ). 
After the best fit was obtained with the experimental data, the theoreti- 
cal dispersion curves were further moved along the f/p-axis to establish 
reasonable limits to the fit. Frcan these limits the probable errors were 
calculated. Values of 9 appear in table II. 


Since it was observed that the dispersion curves gave no evidence 
of the existence of multiple relaxation times, it was assumed that the 
lowest vibrational mode was excited with its characteristic relaxation 
time and that the energy passed freely to the other modes. The relaxa- 
tion time of the lowest mode 9-^ is given by the expression 


9 


1 



( 3*0 


in which is that part of corresponding to the lowest vibra- 

tional mode. Valves of 9^_ were calculated from equation (3*0 and appear 
in table II. 


Collision rates R were calculated from gas kinetic theory, using 
measured or calculated values of gas viscosity, and Z^q was calculated 

for each gas by means of equation (21), using for 9 the relaxation time 
of the lowest mode 9^ calculated from equation (3*1-) . These values of 

Z 10* ’together with corresponding values of R, also appear in table II. 


Experimental Results 

It may be noted in figures 7 through 20 that none of the gases meas- 
ured showed multiple dispersive regions. This indicates that intermodal 
coupling is relatively strong in these molecules, and the assumption made 
by the use of equation ( 3 k) is justified. 

Furthermore, since the use of f/p along the axis of abscissas yields 
the proper form of dispersion curve, the excitation and deexcitation of 
molecular vibrations must be brought about principally by binary collisions. 
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Following the result of Landau and Teller given "by equation (27), 

2/3 1/3 - 

Z]_q has "been plotted against M in figure 21. A theoretical 

curve has "been drawn to fit the data using s = 0.25A and a = 0.5. 

This value of s is about one-eighth the diameter of the colliding mole- 
cules and therefore is not reasonable. A simi l ar curve drawn with 
s = 1A will not fit the data regardless of the value of a used. Fur- 
thermore, the data do not appear to fit a single curve of the type drawn. 
It therefore appears that the theory of Landau and Teller does not fit the 
deexcitation of molecular vibrations in these gases. 

In figure 22 Z 1Q has been plotted against v in accordance with 

equation (29) and the energy-excitation theory. Curves have been drawn 
using the values a = 0.1 and n = 1.8, 1.5, 2.8, and 3.8. The value 
of n to be used apparently depends upon the number of hydrogen atoms 
replaced from the methane molecule. The quantity nhv^ gives the mini- 
mum energy of approach necessary to excite or deexcite the vibrations of 
the molecule. This is presumably the energy necessary for the molecules 
to penetrate deeply into each other’s force fields or to form a compound 
molecule as previously suggested. 

Fogg, Hanks, and Lambert (ref. 12, 1953) have measured the relaxa- 
tion times of some of the halogen-substituted methanes. Their values 
for Z 10 were computed from their relaxation data by making the assump- 
tion that hv^ » kT and hence that Z 10 *=> N0-[_. At values of 

v^_ < 500 centimeter -1 , this is not true and leads to error in Z 10 at 
low frequencies. In accordance with equation (21), the authors' values 

for Z 10 have been multiplied by a factor 1 + e _ ^- v l/ kT . The values 
of Z^o thus corrected are recorded in table HI. 

Since their data were taken at temperatures of about 373 ° K, the 
temperature correction predicted by equation (29) has been used to cor- 
rect Z^o to 300° K; these values are also recorded in table III. The 
corrections thus calculated, however, were invariably larger than the 
differences between these points and the calculated Zj_q curves. The 

values of Z^q for 373° K and also the values predicted for 300° K are 
shown in figure 22. 

Several other authors have measured relaxation times in some of the 
halogen- substituted methanes. Byers (ref. lA, 19^3 ) > Eucken and Aybar 
(ref. 15, 19^0), and Petralia (ref. 16, 1952) have used acoustic inter- 
ferometers to measure relaxation times; Griffith (ref. 17, 1950) and 
Huber and Kantrowitz (ref. 18, 19^7) have made measurements using impact 
nozzles of the Kantrowitz type. Their data are also reported in table III, 
and the corresponding values of Z^q are plotted in figure 22. Where the 
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measurements were made at temperatures more than a few degrees, removed 
from 300 ° K, the data are corrected to this temperature by use of equa- 
tion (29) • 

The fact that use of equation ( 29 ) gives too large temperature cor- 
rections might suggest that the steric factor a is a function of tem- 
perature. It is reasonable to assume that a, being a geometrical 
probability, is dependent upon the speed at which the molecules are 
rotating, which in turn is a function of temperature. The period of a 

-12 

rotation of a typical molecule' is 10 second, which is of the same 
order as the duration of a molecular collision. 


CONCLUSIONS 


Measurements were made of sound dispersion in a number of heavy 
gases. On the basis of, the data reported, the following conclusions 
may be drawn: 

1. Intermodal coupling is so strong in the halogen-substituted 
methanes that once the vibrational mode of lowest frequency is excited 
the energy passes readily into the other modes. 

2. Only binary collisions are important in the excitation of vibra- 
tions in the halogen- substituted methanes within limits of observation. 

3 . The probability, in a given collision, of exciting or deexciting 
molecular vibrations in the halogen-substituted methanes appears to 
depend upon the relative energy of approach of the colliding molecules 
rather than upon their relative velocity as the Landau-Teller theory 
predicts. In particular, the probability of excitation or deexcitation 

is negligible unless the relative energy of approach e exceeds e* = nhv^ 

where e* is the minimum value of e required for vibrational excitation, 
n is characterized by the number of hydrogen atoms replaced, h is 
Planck's constant, and is the lowest vibrational frequency. 

k. The excitation of vibrations in polyatomic organic molecules may 
be a complex process involving the formation of a compound molecule' with 
a very short lifetime. 

5 . The steric factor probably depends upon temperature through its 
dependence upon rotational velocity. 


Iowa State College, 

Ames, Iowa, August 31 , 195^- 
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TABLE I 

MOLECUIAR WEIGHTS AND HEAT CAPACITIES 


Gas 

Purity, 

percent 

M 

v l' 

cm - -'- 

c l> 

eal/mole °K 

^vib* 

cal/mole °K 

CH^Br 

99-4 

94.95 

6ll 

1.03 

2.26 

chcif 2 

97 

86.48 

369 

1.45 

5.99 

CH3CI 

99.5 

50. ^9 

732 

•78 

1.84 

CCI3F 

99.9 

137.38 

244 

4.05 

10.63 

CBrClF 2 

99-9 

165.38 

200 

1.85 

10.24 

CF^ 

— 

88.01 

437 

2.93 

6.91 

C By? 2 

99. 7 

52.03 

532 

1.19 

2.32 

ch 2 cif 

99-9 

68.48 

358 

1.56 

3.46 

CC1F3 

99-0 

104.47 

356 

3.16 

8.16 

CHClgF 

99-0 

102.93 

274 

1.72 

6.76 

CHF3 

— 

70.02 

508 

2.48 

4.78 

CClgFg 

97 

120.92 

260 

1.76 

8.63 

CBr^Q 

99-9 

209.84 

165 

1.88 

10.55 

CBrF* 



148.93 

297 

4.24 

9.84 











TABLE II 


RELAXATION TIMES AND COLLISION LIFETIMES 


OftB 

0, sec 

0^, aeo 

R, sec - ^- 

^O 

V 

ca~^ 

v^V' 5 

CH^Br 

7.5 * 0.3 X 10" 8 

3.4 ± 0.1 x 10 -8 

IO.79 * 0-22 X 10 P 

386 ± 14 

6n 

328 

CHOIFg 

9.5 ± 0.7 

2.3 ± 0.2 

9.98 ± 0.20 

268 t 24 

369 

227 

CHjCl 

20.4 ± 1.1 

8.6 ± 0.5 

10.91 * 0.22 

966 ± 67 

732 

300 

CCljT 

4.0 ± 0.4 

1.3. * 0.2 

11.79 * 0.24 

234 ± 31 

244 

202 

CBxCU^ 

5.0 t o.4 

.9 * ?.l 

9.66 ± 0.48 

121 t 13 

200 

188 

CF4 

73.6 t 0.8 

32.1 ± 0.3 

6.84 * 0.34 

2,460 ± 125 

437 

256 

CE ^ 2 . 

3.0 ± 0.3 

1.5 ± 0.2 

9. 81 t 0.49 

159 1 24 

332 

245 

CHgClF 

1.0 ± 0.1 

„45 * 0.03 

11.20 t O.36 

60 ± 7 

358 

206 

cci? 3 

23.9 i 0.6 

9.3 * 0.2 

11.28 * 0.23 

1,240 t 37 

356 

237 

cHci^r 

2.3 ± 0.4 

.6 ± 0.1 

11.48 ± 0.57 

89 t 15 

274 

196 

CHFj 

47.9 * 1.8 

24.9 * 0.9 

8.34 ± 0.42 

2,260 * l4o 

508 

263 

CCl^?z 

• 7.8 t 0.6 

1.6 ± o.l 

10.31 * 0.21 

212 ± 14 

260 

201 

CBr^Fg 

2.9 ± 0.3 

.50 ± 0.09 

9.34 - 0.47 

68 ± 8 

165 

179 

CBrFi 

19.6 ± 0.1 

8.4 t o.l 

8.05 ± o.4o 

840 t 43 

297 

239 


ro 

vx 
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TABLE in 


BATA CF OTHER INVESTIGATORS 


ro 

P- 


Gas 

(a) 

Investigator 

°K 

6, sec 

^10 

ZiQ at 
T a 300° K 

E9 

Vl 2/ V/ 5 

CH^I 

Fogg, Hanks, and Ianbert (ref. 12) 

373 

4.5 x 10- 8 

86 

176 

533 

542 

CHClj 

Fogg, Hants, and Lambart (ref. 12) 

303 

2.8 

83 

139 

260 

200 

CHgClg 

Fogg, Hanks, and Lambert (ref. 12) 

373 

8.4 

204 

486 


192 

CHjT 

Fogg, Hanks, and Lambert (ref. 12) 

373 

320 

6,070 


1,048 

336 

CCl^ 

Fogg, Hanks, and lambart (ref. 12) 

373 

8.3 

191 

548 

218 

197 

CHjBr 

Fogg, Hanks, and Lambart (ref. 12) 

373 

10 

202 

468 

611 


OHgPg 

Fogg, Hanks, and Lambert (ref. 12) 

373 

5.5 

136 

278 

552 

248 

CHPj 

Fogg, Hanks, and Lambert (ref. 12) 

373 

4a 


3,720 

506 

263 

CHtCI 

Fogg, Hanks, and lambart (ref. 12) 

373 

18 

424 


732 

300 


Griffith (ref. If) 

£95 

20.2 

956 





Petralla (,ref. lb) 

287 

20.2 

93b 




iiiSEMgfc 

Fogg, Hanks, end Lambert (ref. 12) 

373 

66 

1,580 

6,210 

437 

256 


Byers (ref. 1*0 

293 

76 

2,480 




C% 

Euckan and Aybar (ref. 15) 

582 

84 



1,506 

302 

CCl^Pp 

Huber and KantrowitB (ref. 18) 

297 

2 

49 


260 



Griffith (ref. 17) 


8 

218 




CHClTg 

Griffith (raf. 17) 

294 

10 

282 


369 

227 


^H^I, aethyl iodide; CHCI3, irtohlorcmathana; CH0CI2, diC-hiorcnet ham ; CHjE, fluorcraetisne; COl^, tatrachiorcmethane j 
CH^, methane; other ooopoundB identified In text. 
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Figure 1.- Sonic velocity versus frequency for thermal relaxation process. 
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Figure 2.- Interferometer cup shoving crystal mo unting . 



Figure 3*~ Crystal oscillator circuit. 
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Figure 4 .- Crystal impedance -measuring circuit. 


12 h 6AS7 




Figure 6 .- Gas-handling and vacuum equipment. 
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Figure 8.- Velocity dispersion in CHCIF2 
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Figure l4.- Velocity dispersion in CH 2 CIF. 
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Figure 19*- Velocity dispersion in CBr^Fg 



Figure 20.- Velocity dispersion in CBrF- 
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